FEMS Microbes

Antimicrobial resistance (AMR) is an ever-increasing threat to population health. Industrial, environmental and societal factors are increasingly recognized as important contributors to AMR within communities. Here, we investigated the spatial distribution of AMR genes (ARGs) across Alberta, Canada and their association with socio-economic, immigration-related, and agro-industrial characteristics using municipal wastewater-based surveillance. We analyzed monthly wastewater metagenomes collected between March 2022 and March 2023 across eleven municipalities, representing 39% of Alberta's population. Integration with census data enabled multivariate analysis, revealing that municipal resistome profiles were strongly structured along income and immigration-related population gradients. ARGs spanning 14 resistance classes exhibited distinct distributional patterns across income and immigration gradients, including contrasting associations among beta-lactam, aminoglycoside, and macrolide-lincosamide-streptogramin ARGs, consistent with heterogeneous selection pressures across sub-populations. These findings demonstrate the capacity of longitudinal wastewater surveillance to identify persistent population-level resistome patterns and highlight the importance of incorporating sociodemographic context into AMR surveillance and mitigation strategies.

Objectives: To compare the effectiveness of povidone-iodine ear wash compared to no ear wash and oral cotrimoxazole compared to placebo given in addition to standard topical antibiotic treatment (ciprofloxacin drops) for chronic suppurative otitis media (CSOM) in Australian Aboriginal children. Methods: A randomised, parallel, 2 x 2 factorial design, assessor-blinded clinical trial in the remote Northern Territory of Australia. Aboriginal children with confirmed CSOM were eligible to be randomised into four treatment groups, allowing two primary treatment comparisons in a 2-in-1 trial approach. Participants received standard treatment (twice daily cleaning and topical ciprofloxacin drops) plus: i) either 16 weeks of pre-treatment povidone-iodine ear wash or no povidone-iodine ear wash; and ii) either 16 weeks of oral cotrimoxazole or placebo. Central randomisation with allocation concealment and triple-blinding of the oral antibiotic treatment arms was used. The relative risk (RR) and risk difference (RD) were estimated after adjustment for age, community, and the other intervention. The primary outcome was the proportion of children with any otorrhoea (clinical failure) after 16 weeks of treatment. Secondary outcomes included size of tympanic membrane (TM) perforation and amount of discharge, time to cessation of discharge, proportion of children with respiratory and other pathogens in ear discharge (at baseline and 16 weeks) and hearing levels (at 12 months). Findings: 280 children with CSOM were randomised and 270 had their primary outcome assessed. Clinical failure (presence of any ear discharge) after 16 weeks of treatment was 66/134 (49%) in the povidone-iodine group versus 69/136 (51%) in the no povidone-iodine group (RD= -1% (-12,11), p= 0.93) and 56/134 (42%) in the cotrimoxazole group versus 79/136 (58%) in the placebo group (RD=-16% (-28,-4), p=0.007). The amount of discharge, TM perforation size, the level of hearing impairment, and serious adverse events were not significantly different in both treatment comparisons. Anaerobic growth (24%), Pseudomonas aeruginosa (21%) and Haemophilus influenzae (17%) were the most common pathogens found in the ear discharge before treatment. Fungi or yeast (24%), Staphylococcus aureus (15%), and anaerobic growth (10%) were the common pathogens after 16 weeks of treatment, with no significant differences between groups. At 12 months post-randomisation, 55-60% of children had at least one discharging ear and there was no difference between treatment groups. Interpretation: Povidone-iodine ear washes did not contribute to better ear outcomes in this study. Cotrimoxazole for 16 weeks resulted in more children with clinical improvement to dry ears. Oral cotrimoxazole may play a role in reducing the burden of CSOM in populations with high rates of persistent disease.

Background Coccidioidomycosis is an emerging fungal disease across the arid Americas and a frequent cause of community-acquired pneumonia. Understanding where Coccidioides populations originate, how they move across space, and whether they are expanding is important for interpreting changing patterns of Valley fever and anticipating future infection risk. Methods We prospectively collected and whole-genome sequenced 186 Coccidioides-positive clinical isolates submitted to a national diagnostic laboratory, and included 126 previously sequenced genomes. We applied genomic clustering, time-calibrated phylogenetic reconstruction, ancestral area reconstruction, mating-type assignment, and demographic inference to identify major populations, infer dispersal patterns, assess evidence for recombination and clonality, and reconstruct historical population dynamics. Findings We analyzed 312 genomes (139 C. immitis; 173 C. posadasii) and identified three major genetic populations within each species. C. immitis included two California-centered populations and one Pacific Northwest population, whereas C. posadasii included two Arizona-centered populations and one Texas-centered population. The most recent common ancestor was estimated at approximately 127,000 years for C. immitis and 234,000 years for C. posadasii. Most populations were not fully monophyletic, consistent with retained ancestral variation and/or ongoing gene flow. Inferred dispersal was largely asymmetric, with most movement originating from California in C. immitis and from Arizona and Texas in C. posadasii. Most populations contained both mating types, but one C. immitis population and a Brazilian subgroup of C. posadasii were clonal. All populations showed recent demographic expansion. Interpretation The evolutionary history of Coccidioides is characterized by strong geographic structure, ongoing gene flow, and recent demographic expansion. These processes are likely to influence future patterns of Valley fever endemicity and supports the use of genomic surveillance to detect shifts in disease risk as environmental conditions change.

Exposure to higher levels of ambient air pollution during pregnancy has been linked to multiple adverse pregnancy outcomes. However, studies on acute exposures and reduced gestation length have reported inconsistent findings. This project aims to examine the acute association between ambient air pollution and preterm (28-36 gestational weeks) or early-term (37-38 gestational weeks) births. Daily concentrations of 12 air pollutants, based on bias-corrected numerical model outputs, were linked to vital records of singleton live preterm and early-term births from 2005-2017 in California, Florida, Georgia, Kansas, Nevada, New Jersey, North Carolina (2005-2015), and Oregon. Under a time-stratified case-crossover design, odds ratios (OR) were estimated via conditional logistic regression with adjustment for risks among ongoing pregnancies, meteorology, time trends and federal holidays. We estimated cumulative associations up to a 6-day lag using distributed lag models. Risk estimates per interquartile range (IQR) increase in exposure were pooled across states using inverse-variance weighting. Our study included 1,085,162 preterm and 3,901,185 early-term births. We observed positive associations between 0-2 day cumulate exposure to several air pollutants and early-term births, including NO2 (OR: 1.0023, 95% CI:1.0010, 1.0037 per 7.1 g/m3 increase), PM2.5 (OR=1.0022, 95% CI: 1.006, 1.0038 per 4.6 g/m3 increase), PM2.5 organic carbon (OR= 1.0026, 95% CI: 1.0013, 1.0039 per 1.7 g/m3 increase) and PM2.5 elemental carbon (OR=1.0025, 95% CI: 1.0014, 1.0035 per 0.26 g/m3 increase). Associations with preterm birth were mostly null. In conclusion, we found positive associations between short-term air pollution exposure, including PM and major PM2.5 components, and risks of early-term birth.

Background: Early-onset colorectal cancer (EOCRC) is increasing globally, particularly among Hispanic/Latino (H/L) populations, yet the contribution of tumor-colonizing microbiota to age-associated colorectal cancer (CRC) biology remains poorly understood. Most microbiome studies have focused on fecal communities or non-Hispanic populations, leaving the intratumoral microbial landscape of H/L patients largely unexplored. Methods: We performed an exploratory characterization of tumor-colonizing microbiota using whole-exome sequencing (WES) data from four primary colorectal tumors obtained from H/L patients treated at City of Hope, including two EOCRC (<50 years) and two late-onset colorectal cancer (LOCRC; [&ge;]50 years) cases. Following removal of host-derived sequences, microbial taxonomic profiling was conducted at the family, genus, and species levels, and microbial metabolic pathways were inferred. Clinical and pathological data were integrated to evaluate age-associated differences in microbial composition and predicted function. Results: Family-, genus-, and species-level analyses consistently demonstrated greater microbial diversity in LOCRC than EOCRC. LOCRC contained more than twice the number of unique bacterial families, nearly three times as many unique genera, and more than twice as many unique bacterial species. A conserved core microbiota, including Fusobacteriaceae, Prevotellaceae, Fusobacterium, and Prevotella, was identified across both age groups, whereas LOCRC was enriched in CRC-associated taxa including Fusobacterium nucleatum, Bacteroides fragilis, Parvimonas micra, Porphyromonas asaccharolytica, and Dialister pneumosintes. Species-level analyses revealed only a single shared bacterial species between EOCRC and LOCRC, indicating progressive microbial divergence with increasing taxonomic resolution. In contrast, functional profiling identified 11 predicted microbial metabolic pathways, of which nine were shared between age groups, two were unique to EOCRC, and none were exclusive to LOCRC. Core metabolic pathways involved in energy metabolism, amino acid biosynthesis, phospholipid metabolism, and central carbon metabolism exhibited comparable abundance across both groups, demonstrating substantial functional conservation despite pronounced taxonomic differences. Conclusions: Tumor-colonizing microbiota differ markedly between EOCRC and LOCRC in H/L patients, with late-onset tumors exhibiting substantially greater microbial richness and taxonomic complexity. Despite these compositional differences, microbial metabolic functions remain largely conserved, supporting the concept of functional redundancy within the colorectal tumor microenvironment (TME). Although exploratory, this proof-of-concept study provides one of the first characterizations of intratumoral microbiota in H/L EOCRC and establishes a foundation for larger multi-omics investigations aimed at identifying microbiome-based biomarkers and therapeutic targets for precision oncology.

Blood banks represent highly controlled healthcare environments where microbiological surveillance has traditionally focused on blood products rather than environmental microbial reservoirs. Despite their critical role in transfusion safety, the ecology of surface-associated microorganisms and the persistence traits that enable their long-term survival remain poorly understood. Here, we combined scanning electron microscopy, culture-based microbiology, phenotypic characterization, MALDI-TOF mass spectrometry, and whole-genome sequencing to investigate whether surfaces within a public blood bank facility constitute reservoirs of environmentally derived bacteria with enhanced persistence potential. Samples collected from a public blood bank in Tucuman, Argentina yielded 37 culturable bacterial isolates, predominantly Gram-positive environmental taxa together with a limited number of opportunistic Gram-negative species. More than 30% of the isolates exhibited multidrug resistance, while several strains displayed strong biofilm formation, amyloid-like fiber production, motility, and hemolytic activity, indicating multiple phenotypic strategies associated with long-term surface persistence. Whole-genome sequencing of six representative isolates confirmed species identity, identified genes related to antimicrobial resistance, adhesion, biofilm formation, stress adaptation, and cytotoxicity, and revealed frequent genotype-phenotype discordance, highlighting the importance of integrating genomic and phenotypic analyses. Notably, one isolate exhibited less than 92% average nucleotide identity with publicly available genomes, suggesting the presence of a previously undescribed environmental species. Thus, blood bank surfaces function as selective ecological niches favoring bacteria with persistence-associated traits rather than simply reflecting contamination from blood products. These microorganisms may constitute latent biosafety hazards if environmental barriers fail, particularly in facilities handling biological materials intended for vulnerable patients. Our results support the incorporation of integrated bioimaging, phenotypic characterization, and genome-resolved environmental surveillance into infection prevention strategies and transfusion biosafety programs within a One Health framework.

Diarrheal disease remains a leading cause of morbidity and mortality among children under five globally, with poor sanitation and hygiene accounting for over 88% of diarrhea and malnutrition burden. In Kenya, diarrhea ranks third in under-five mortality, particularly affecting arid and semi-arid regions. School-Led Total Sanitation (SLTS), adapted from Community-Led Total Sanitation (CLTS), uses pupils as change agents for household hygiene knowledge transfer. However, SLTS effectiveness in addressing diarrhea and malnutrition has not been evaluated in Kenya. This study assessed SLTS effects on diarrheal disease and nutritional outcomes among children aged 5-59 months in Baringo County. A pre- and post-test quasi-experimental design with nonequivalent control groups was employed in Mogotio (intervention) and Baringo South (control) sub-counties. Using multistage sampling, 440 children aged 6-59 months were enrolled. The six-month SLTS intervention included capacity building, school health club formation, triggering activities using Participatory Rural Appraisal tools, Information, Education, and Communication materials distribution, continuous sensitization, and household monitoring. Data were collected at baseline and three months post-intervention using electronic questionnaires and anthropometric measurements. Nutritional status was assessed using WHO Anthro software z-scores for length/height-for-weight (HWZ), and weight-for-age (WAZ) to determine wasting, and underweight prevalence. Chi-square analysis assessed intervention-control differences. Baseline and endline socio-demographic characteristics were comparable between groups. At endline, no significant nutritional outcome differences were observed: wasting prevalence was at 15.0% versus 16.4% ({chi}{superscript 2}=0.155, df=1, p=0.694) while underweight was 12.3% versus 13.6% ({chi}{superscript 2}=0.181, df=1, p=0.670). However, diarrheal disease prevalence significantly reduced in intervention versus control groups: 5.9% versus 13.2% ({chi}{superscript 2}=6.738, df=1, p=0.009), representing a 53% reduction. SLTS intervention showed no significant effect on nutritional outcomes but demonstrated a significant reduction in diarrheal disease among children aged 6-59 months. These findings provide strong evidence for integrating school-based sanitation and hygiene approaches into broader public health strategies addressing diarrheal diseases.

Background and Aims The optimal treatment for first episodes and first recurrences of Clostridioides difficile infections (CDI) is unknown and there is emerging evidence for pulse and taper (P-T) regimens. Therefore, we sought to estimate the relative efficacy of treatment options. Methods MEDLINE and CENTRAL were searched from database inception to May 21, 2025 and unpublished conference abstracts were searched from recent infectious disease conferences. RCTs on the treatment of first episodes or first recurrences of CDI comparing fixed-dose or P-T regimens of fidaxomicin or vancomycin were included. The primary and secondary outcomes were 40- and 56-day CDI recurrence, respectively. A random-effects network meta-analysis on the risk ratio (RR) scale was conducted using a standard regimen (10-14 days) of vancomycin as the comparator. Treatments were ranked using the surface under the cumulative ranking curve (SUCRA). Results 8 RCTs were included comprising a total of 2181 patients. For 40-day recurrence, fidaxomicin P-T had the highest probability of ranking best (RR=0.10, 95%Confidence Interval [95%CI]=0.10-0.49, SUCRA=1.00), followed by vancomycin P-T (RR=0.49, 95%CI=0.32-0.76, SUCRA=0.61), fixed-dose fidaxomicin (RR=0.61, 95%CI=0.49-0.76, SUCRA=0.39), and, finally, fixed-dose of vancomycin (SUCRA=0.00). The treatments ranked in the same order for 56-day recurrence, though only 3 RCTs reported on this timepoint. Conclusion Vancomycin P-T, fidaxomicin P-T, and fixed-dose fidaxomicin were all superior to a fixed-dose vancomycin. Head-to-head comparative effectiveness RCTs are needed to quantify their relative effect sizes of and impact on long-term prevention of recurrent CDI.

The World Health Organization (WHO) validated Sri Lanka had eliminated lymphatic filariasis as a public health problem in 2016, the second country in Southeast Asia to attain this status. However, post-validation surveillance has identified sporadic cases of brugian filariasis. The reemergence of Brugia malayi infections in Sri Lanka warrants urgent investigations. Recent studies have shown that the parasite responsible for the reemergence is a novel zoonotic Brugia sp. maintained among dogs that is closely related but distinct to the human-infecting B. malayi species. The current study employed morphological and morphometric assessments, revealing that this novel zoonotic Brugia sp. is within the B. malayi morphological range. Molecular characterization of three genomic regions, the nuclear genomic region SLXI, the non-coding region HhaI, and the mitochondrial genomic region COXI confirmed it as a genetic variant more closely related to B. malayi than to B. pahangi. Phylogenetic analysis further indicated it as a distinct genomic variant, closely related to a B. malayi-like parasite reported from India. Notably, that same parasite was identified in infected humans, animals, and potential vector mosquitoes. This, together with the detection of both human and animal blood within the same brugian infective mosquitoes, and delineating the canine origin of the parasites in human infections, provides compelling evidence supporting zoonotic transmission of this parasite. To our knowledge, this is the first report demonstrating the presence of the same brugian parasite in humans, domestic animals, and potentially infective mosquitoes in Sri Lanka, supported by multi-genomic evidence. The recent identification of multiple potential mosquito vector species suggests that this parasite may have undergone adaptive changes, facilitating its ability to overcome the species barrier. These findings substantiate the long-held hypothesis of zoonotic transmission of the reemerged brugian parasite, highlighting significant implications for ongoing surveillance and control strategies.

Background: Monoclonal antibodies have emerged as a prophylactic strategy to prevent symptomatic SARS-CoV-2 infection in immunocompromised individuals. However, the evolutionary and clinical implications of breakthrough infections under this regime remain unclear. Methods: A male in their 80s with a haematological/oncological diagnosis received a 2000 mg intravenous infusion of sotrovimab in March 2023 and was diagnosed with COVID-19 by RT-qPCR from a nasopharyngeal swab in August 2023. Weekly samples (n=24) were collected through February 2024 (171 days). All samples underwent whole-genome sequencing, with select mutations subjected to functional assessment. Findings: Sequencing identified the GE.1 lineage at all timepoints. An intra-host recombination event in ORF1ab (positions 8942-12458) was detected prior to 23 weeks post-detection, followed by a 14-fold increase in viral load (7.42e+06 to 1.00e+08 RNA copies/mL) and a marked shift in the viral population. E340D, a sotrovimab resistance mutation, was detected at low abundance (46%) within the first week post-infection, fluctuated over time, and was nearly fixed by week 15 (107 days) post-detection. We assessed five spike mutations - V36M, S98F, and V213G in the N-terminal domain, Y505P in the receptor-binding domain, and P681Q near the S1/S2 cleavage site - and additionally evaluated the impact of E340D. V36M conferred the highest infectivity across all cell lines, with the most significant effect in low-TMPRSS2 cells. While all mutations showed enhanced infectivity with the addition of E340D, the effect was most pronounced in mutations with lower baseline infectivity. The addition of E340D significantly decreased relative neutralizing titres for V36M, S98F, and V213G, enabling escape from neutralizing antibodies in XBB-responsive individuals, illustrating an enhanced phenotypic advantage. Patient neutralizing activity was absent pre-sotrovimab, and sotrovimab-induced neutralization was further compromised by selection of E340D. Interpretation: Sotrovimab pre-exposure prophylaxis in an immunocompromised patient did not prevent SARS-CoV-2 infection, and selected for resistant mutation E340D, with unexpected fitness consequences across non-receptor binding domain spike regions.

We present FootNet, a 453-image multi-view smartphone foot dataset for binary foot segmentation, with expertannotated masks across six anatomical views (dorsal, medial, and plantar, both left and right). We benchmark four segmentation models under a controlled protocol: U-Net with a MobileNetV2 encoder achieves the best performance (IoU 0.9268, Dice 0.9608, 95 % CI [0.9209, 0.9320]); DeepLabV3 with MobileNetV3-Large scores IoU 0.8984 (Dice 0.9449); UNet++ with MobileNetV2 scores IoU 0.8913 (Dice 0.9391); and SAM ViT-B with oracle boundingbox prompt scores IoU 0.9219 on the matched 191-image subset. Bonferroni-corrected Wilcoxon signed-rank tests (k = 6 comparisons) show U-Net significantly outperforms DeepLab (p < 0.001, r = 0.638) and SAM ViT-B with oracle boundingbox (p = 0.005, r = 0.202); UNet++ does not significantly differ from DeepLab (p = 0.062). Connected-component postprocessing yields negligible benefit (mean {triangleup}IoU = +0.0003, 12 of 453 images improved). The extended dataset is available upon request

Purpose: Nonresponse to routinely collected patient-reported outcome measures (PROMs) threatens the representativeness of aggregated data. We characterized patient-, provider-, and clinic-level factors associated with PROMIS Global-10 nonresponse in routine radiation oncology care. Methods: In this retrospective cohort study, all adults seen at five Mass General Brigham radiation oncology clinics over one year were included. The primary outcome was patient-level nonresponse, defined as never completing the portal-administered Global-10 versus completing it at least once. Using iterative mixed-effects logistic regression, we modeled patient-, provider-, and clinic-level factors. Results: Among 12,214 patients, 71 providers, and five clinics, patient- and appointment-level response rates were 35.4% and 10.9%, with patient-level response ranging nearly fivefold across clinics (12.8% to 66.2%). In Model 1, male sex, lower education, not working, and recent surgery had higher odds of nonresponse, and longer time since diagnosis lower odds. After provider- and clinic-level factors were added, patient sex, education, and employment became nonsignificant, whereas recent surgery (adjusted odds ratio [aOR] 1.97) and longer time since diagnosis (aOR 0.46 for >12 months) persisted. A provider's historical collection rate was protective but attenuated at the clinic level. There, a later program launch (aOR 0.29) and higher historical collection rate (aOR 0.79) correlated with lower nonresponse, whereas academic versus community setting did not. Conclusions: Nonresponse to routinely collected PROMs is a multilevel phenomenon driven substantially by clinic-level implementation factors, not patient characteristics alone. Because response rate is only a proxy for representativeness, PROMs programs and PRO-based performance measures should prioritize representative collection over volume.

Controlled human infection models (CHIMs) are increasingly used in infectious disease research to study pathogen dynamics and evaluate interventions under controlled conditions. However, these studies are resource-intensive and involve ethical and safety constraints, making efficient study design critical. Dose-finding is a key early component in CHIMs, where the aim is to identify a challenge dose that achieves a target infection probability. Traditional rule-based designs are commonly used but can be inefficient, motivating the use of model-based adaptive approaches such as the Bayesian Continual Reassessment Method (CRM). Although CRM has been extensively studied and widely adopted in Phase I oncology trials for identifying the maximum tolerated dose of therapeutics, its application in CHIM settings remains limited, particularly when the endpoint of interest is infection. This tutorial provides step-by-step guidance for implementing a Bayesian CRM in dose-finding CHIMs, using an oropharyngeal Neisseria gonorrhoeae challenge as a motivating case study. The framework outlines key design components, including dose-grid specification, dose-response model, prior elicitation, Bayesian updating, decision rules, and stopping criteria, with particular emphasis on a clinically interpretable parameterisation. Trial operating characteristics are evaluated through simulation studies under multiple dose-response scenarios and prior-predictive analyses, and compared with a commonly used '3+3' type rule-based design. This work highlights the advantages of Bayesian model-based designs for dose-finding in CHIMs over classic rule-based designs and provides a structured, reproducible framework for implementing CRM, supporting their application in future CHIM studies.

Background The 2022-2023 global mpox outbreak highlighted the critical need for robust genomic surveillance capabilities to track mpox virus (MPXV) evolution and transmission dynamics. Methods Building upon our established SARS-CoV-2 sequencing infrastructure, we implemented a Molecular Loop probe-based long-read sequencing approach using Pacific Biosciences Sequel II technology for comprehensive MPXV genomic surveillance across the United States (US). From August 2024 to June 2025, we generated 326 high-quality whole genome sequences from residual mpox-positive clinical specimens collected by Labcorp across all 10 US Department of Health and Human Services regions. Results Our analysis identified two samples containing clade Ib MPXV in January and June 2025 and captured shifting trends in clade IIb diversity, with 13 distinct lineages observed. We also identified multiple instances of large (~1.6-17.6kb) deletions proximal to the inverted terminal repeats in clade IIb genomes. APOBEC3 mutation analysis indicated substantial evidence of human-to-human transmission among both clades. Further, we observed significantly higher APOBEC3-associated SNPs per kilobase (P<0.001) in clade IIb genomic variable regions relative to their central conserved region. Our assay exhibited strong reproducibility across biological replicates from individual patients and accuracy was confirmed via parallel sequencing of select specimens by US Centers for Disease Control and Prevention (CDC) using metagenomic sequencing. We also demonstrated via custom simulation that our assay discriminates all known MPXV clades and lineages, including those we have not observed in the US. Conclusions Our integrated nationwide surveillance system facilitates real-time genomic tracking of outbreak evolution, with demonstrated capacity across SARS-CoV-2 and MPXV, positioning this platform for rapid deployment during future pathogen emergence.

Stochastic epidemic models are a cornerstone of infectious disease epidemiology and are often used to study intervention scenarios. However, large run-to-run variability can make intervention effects difficult to estimate precisely. We revisit the epidemic Sellke construction, which assigns each individual an infection threshold for the cumulative infection hazard such that, conditional on the thresholds, the epidemic trajectory becomes deterministic. This enables coupling of simulations with and without an intervention, yielding low-variance effect estimates even when outcomes such as final size or peak incidence vary widely between runs. We develop an exact, event-driven implementation that maintains infection and recovery events in priority queues. Cumulative infection-hazard updates require O(log N) time per event, yielding overall complexity O(Elog N) for E events in a population of size N. The implementation achieves computational performance comparable to the classical Gillespie algorithm while naturally accommodating non-Markovian infectious periods and complex infectiousness profiles. We illustrate the approach using distance-dependent spread of avian influenza between poultry farms in the Netherlands and a multilayer population with households, schools, and workplaces. In both examples, coupling enables efficient within-run comparisons of intervention scenarios across stochastic realisations.

Background: Opioid use disorder (OUD) is heritable, yet most genome-wide association studies (GWAS) have focused on European populations, leaving the genetic architecture of OUD in non-European populations underexplored. Methods: We conducted GWAS of OUD across three ancestries using electronic health records and genomic data from 52,357 All of Us Research Program participants (8,912 cases; 43,445 matched opioid-exposed controls; 48.5% female). Participants were stratified into European (EUR), African (AFR), and Admixed American (AMR) ancestry groups for logistic regression GWAS, with independent replication in the Million Veteran Program. We then applied the deep-learning model AlphaGenome to predict the tissue-specific transcriptomic and splicing consequences of top risk variants across 13 reward-pathway brain regions. Results: We identified and replicated a novel DDX6 risk locus, alongside established OPRM1 and FURIN signals. AlphaGenome predicted the DDX6 regulatory allele downregulates the stress-resistance gene FOXR1 in the nucleus accumbens, while the protective OPRM1 variant (rs1799971) upregulates OPRM1 expression across reward networks. Other signals of interest included IL6R and SHISA9 (EUR); GHR (AFR); and ASTN2 (AMR). Conclusions: This study identifies DDX6 as a novel OUD risk locus, replicates associations with OPRM1 and FURIN, and highlights biologically plausible ancestry-specific signals in AFR and AMR populations. We also replicated top variants in an independent population. Finally, integrating GWAS with deep-learning annotations provides specific, localized biological hypotheses to guide future experimental validation and targeted therapeutics.

Background: International Classification of Diseases (ICD) codes are often used in epidemiological studies to track disease rates over time. Objective: This evaluation of ICD-code-based algorithms for electronic medical record (EMR) studies of Alzheimers disease (AD) and related dementias (ADRD) examines the impact of incorporating Centers for Medicare and Medicaid (CMS) data as an additional source of diagnostic and treatment information in Department of Veterans Affairs (VA) EMR studies. Methods: We performed a chart review of 100 VA Million Veteran Program (MVP) participants to evaluate algorithm performance. We also assessed genetic associations across algorithms in a large MVP cohort (n=396k). Results: Adding CMS data increased the number of detected cases, sensitivity, and positive predictive value, but decreased specificity and negative predictive value. Genetic analyses showed that broader (ADRD/dementia) algorithms with just VA data performed similarly to narrow (AD-focused) algorithms incorporating both VA and CMS ICD codes. Additionally, narrow AD algorithms based solely on VA data yielded the highest ORs, indicating the largest proportion of late-onset AD cases. Conclusions: We recommend using a broad (ADRD) algorithm without CMS or medication data, particularly for epidemiological studies or a strict AD algorithm including CMS and medication cases for genetic discovery of late-onset AD associations in VA EMR, and a strict AD algorithm without CMS data for applications focused solely on AD and sensitive to misspecification. Careful evaluation of algorithm performance is warranted in different EMR systems, as ICD coding practices vary by institution, as demonstrated by this comparison of VA EMR and CMS data.

Estimating dengue force of infection (FOI) is essential for understanding transmission dynamics and targeting intervention programmes, yet surveillance data in endemic settings required for estimations are often incomplete, with varying formats. We developed a Bayesian hierarchical catalytic model that jointly fits age-stratified case data, aggregate case data, and seroprevalence surveys within a single framework, incorporating external covariates to improve parameter identifiability. Synthetic validation showed that covariates alone recovered accurate FOI point estimates even when most districts contributed only aggregate data, but did so with poorly calibrated uncertainty; anchoring the model with a single seroprevalence survey was necessary to bring credible interval coverage close to nominal. Applied to 128 districts across Java and Bali, Indonesia (2016-2024), the model revealed substantial spatial heterogeneity in FOI and reporting rates. Many districts in Java exceeded the WHO-suggested seroprevalence threshold for vaccine introduction, yet were classified as low-priority when using reported incidence as prioritisation criterion, particularly in areas with weak surveillance. Model-based seroprevalence estimation, integrating multiple data sources, offers a more consistent basis for identifying high-priority districts for vaccine introduction, and is less susceptible to surveillance bias than reported incidence.

Background: Children born prematurely (before 37 weeks) or admitted to the neonatal unit (NNU) are at increased risk of adverse long-term physical health outcomes. It is also recognised that there is an association with later academic performance and special educational needs, however it is not clear whether these broad risk factors could be used as stand-alone heuristics to identify children who may benefit from additional support in educational settings. We aimed to examine the associations between neonatal unit (NNU) admission and educational attainment in mid-childhood. Methods and Findings: Pregnancy data from a prospective birth cohort (Pregnancy Outcome Prediction Study, Cambridge, United Kingdom, 2008-2012) were linked to national educational outcomes (Department for Education, United Kingdom). Multivariable regression models adjusted for maternal, child, and socioeconomic factors were used to evaluate associations between (i) all NNU admissions, (ii) at term NNU admissions >48 hours, (iii) preterm birth without ongoing physical health needs, and educational outcomes at ages 5-11 years. Children who required any NNU care were more likely not to meet expected educational standards across multiple ages and domains in early and mid-childhood: age 5 early year foundation (aOR 1.64, 95% CI 1.19-2.27, p=0.003), phonics at age 6 (aOR 2.43, 95% CI 1.72-3.57, p<0.001), and at age 7 (here assessments were divided into multiple domains): reading (aOR 1.67, 95% CI 1.18-2.38, p=0.004), writing (aOR 1.72, 95% CI 1.25-2.38, p<0.001), mathematics (aOR 1.56, 95% CI 1.09-2.22, p=0.020), and science (aOR 1.85, 95% CI 1.22-2.78, p=0.003). Similar patterns were observed among both at term-born infants who stayed >48hrs in NNU (phonics assessment at age 6 aOR 2.26, 95% CI 1.51-3.36, p<0.001) and in children born preterm without long-term physical health sequelae (phonics assessment at age 6 aOR 3.07, 95% CI 1.96-4.81, p<0.001). These associations were robust to adjustment for demographic, perinatal, and socio-economic factors. By age 11, differences in academic attainment were attenuated and no longer clearly distinguishable across all exposure groups. However, there was an increased likelihood of special educational needs (SEN) at age 11 associated with any NNU admission (aOR 1.78, 95% CI 1.15-2.73, p=0.009), at term NNU admission for >48hrs (aOR 1.88, 95% CI 1.19-3.00, p=0.007), and children born preterm without long-term physical health sequelae (aOR 1.50, 95% CI 1.00-2.25, p=0.049). Predictive performance of any NNU admission for SEN at age 11 was moderate (AUC 0.70, 95% CI: 1.14-2.65, p=0.010), with balanced sensitivity and specificity and high negative predictive value. Conclusions: NNU admission, for both term and preterm infants, is associated with poorer educational outcomes and an increased likelihood of special educational needs in mid-childhood.

Background Post-infection conditions (PICs), such as Long Covid, are associated with heterogeneous, fluctuating symptoms that profoundly affect daily functioning. Despite moderate-certainty evidence from the NIHR-funded LISTEN trial (COV-LT2-0009) that personalised self management support improves outcomes and may reduce societal and economic impacts of Long Covid, many people living with PICs still receive condition-specific services, generic advice, or stand-alone digital tools that do not address their complex needs. Aim To map care approaches in general practice and synthesise UK evidence for PIC management. Design and setting Scoping review and online survey. Method A two-phase study was conducted: (1) a scoping review of UK evidence on PIC management in general practice; and (2) a supplementary online survey of practitioners working in UK general practice to provide contextual insights. Results The scoping review identified 32 studies focused on Long Covid. One study included a comparator group (ME/CFS). Study populations were predominantly white ethnicity and female. Evidence for non-Covid PICs in UK general practice was largely absent. The supplementary survey (n=46) provided preliminary practice-level insights. Healthcare practitioners reported varied PIC presentations, diagnostic uncertainty, limited referral pathways, inequitable access, and low confidence in managing PICs. Conclusion Evidence informing PIC management in UK general practice remains predominantly Long Covid-focused and may not reflect the range of PICs encountered in practice. While survey findings are preliminary and require confirmation in larger samples, they highlight uncertainty around PIC management. Further research is needed to evaluate whether existing Long Covid pathways should be expanded or complemented by broader PIC models. Keywords general practice; Long Covid; self-management; post-viral syndromes